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Abstract 

The internal structure of a charged spherical black hole is still a topic of debate. In a 
nonrotating but aspherical gravitational collapse to form a spherical charged black hole, 
the backscattered gravitational wave tails enter the black hole and are blueshifted at the 
Cauchy horizon. This has a catastrophic effect if combined with an outflux crossing the 
Cauchy horizon: a singularity develops at the Cauchy horizon and the effective mass 
inflates. Recently a numerical study of a massless scalar field in the Reissner-Nord- 
strom background suggested that a spacelike singularity may form before the Cauchy 
horizon forms. We will show that there exists an approximate analytic solution of the 
scalar field equations which allows the mass inflation singularity at the Cauchy horizon 
to exist. In particular, we see no evidence that the Cauchy horizon is preceded by a 
spacelike singularity. 
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1 Introduction 

The final state of a star that undergoes gravitational 
collapse into a black hole is described by the unique- 
ness theorems of general relativity M. At late times, 
after the star's irregularities have been radiated away 
the external geometry is described by the station- 
ary Kerr-Newman solution. But the structure of the 
realistic black hole's interior has not been definitively 
determined and is still the subject of much debate. 

The unproven, yet plausible strong cosmic censor- 
ship principle leads one to suspect that the singularity 
in a physical black hole ought to be spacelike, and de- 
scribed by the general mixmaster type solution ||] . But 
the Kerr-Newman singularity at r — is timelike. How 
could the inclusion of the details of collapse alter this 
picture so drastically? Penrose has described the basic 
mechanism B . In the Kerr-Newman solution (and the 
Reissner- Nordstrom solution), the timelike singularity 
is preceded by a Cauchy horizon. The Cauchy horizon 
is the boundary of predictability, making the existence 
of a timelikc singularity in a region beyond it physi- 
cally irrelevant. The Cauchy horizon is also a surface of 
infinite blucshift. Any freefalling observer crossing the 
Cauchy horizon will measure infalling radiation to have 
infinite energy density. When only ingoing radiation is 
present a weak nonscalar singularity forms, which is 
classified as a whimper singulartity || . Whimper sin- 
gularities are known to be unstable to perturbations 
which can transform them into stronger singularities 
1- 

This scenario can be examined more closely by as- 
suming spherical symmetry. When an outflow cross- 
ing the Cauchy horizon is included in the analysis of a 
charged spherical black hole |7|, ^ , the mass function is 
found to diverge exponentially. In spherical symmetry, 
the only non-vanishing component of the Weyl tensor, 
the invariant ^2, is proportional to the mass function. 
The net result of including an outflux is that a scalar 
curvature singularity forms along the Cauchy horizon. 
This phenomenon has been dubbed "mass inflation" 
H.The Reissner-Nordstrom black hole has the same 
causal structure as the Kerr-Newman black hole, so we 
expect that results for spherical symmetry should be 
qualitatively the same for the more general black hole. 
Preliminary calculations suggest that the general pic- 
ture derived for spherical symmetry does not change 
dramatically in a non-spherical black hole ||. 

The important issue, which will be addressed in this 
paper, is to ask how generic is the mass inflation pic- 
ture. General stability arguments JlO[ and the numeri- 
cal study of a scalar field |ll[ suggest that the Cauchy 
horizon is preceded by a spacelike r = singularity. 
Hence evolution of the interior should end before the 



Cauchy horizon forms. But while these studies are sug- 
gestive, a conclusive analysis of the internal structure 
is still lacking even for spherical holes. 




Fig 1. Penrose conformal diagram of a col- 
lapsing star. Note that the point H is not part 
of the manifold, but a singular point of this 
mapping. 

Our main objective is to examine the "corner" re- 
gion, ("point" H in Fig. 1), bounded by the Cauchy 
and event horizons, and to decide whether it is possi- 
ble for an early portion of the Cauchy horizon to sur- 
vive. We have found, using analytic approximations 
and physical arguments, that there is a general solu- 
tion of the scalar field equations which is consistent 
with the mass inflation picture [ fL2[ . 

The description of the black hole interior is simplified 
by causality: behind the event horizon the coordinate 
r is timelike, so a descent into a black hole is an evo- 
lution in time. The solution down to any particular 
radius is only influenced by the initial data at larger 
radii. Most importantly, our lack of understanding of 
the true quantum description of gravity (applicable in 
the innermost layers) will not affect the description of 
the outer classical and semi-classical layers. The clas- 
sical evolution problem is thus on firm ground. 

In an aspherical collapse with zero angular momen- 
tum, the asphericities must be radiated away in order 
to form a static spherically symmetric black hole, as 
the no hair theorems demand. Price |2J] has shown 
that the star settles into an asymptotically Reissner- 
Nordstrom state with perturbations that die off with 
advanced time as an inverse power law. These per- 
turbations are scattered by the black hole's external 
potential barrier so that part of the radiation falls into 



the black hole. Near the event horizon, the backscat- 
tered energy flux takes the form v~ p where v — > oo 
corresponds to X + and p = Al + 4, where I is the spher- 
ical harmonic multipole of the perturbation |13| . This 
result was recently verified by numerical studies of the 
coupled Einstein-scalar field equations [Q. We will 
use the Price tail v~ p as the initial data on the event 
horizon, the natural choice for the initial data surface 
for the black hole interior evolutionary problem. 

Our paper is organised as follows. In §2 and §3 we 
will review the mass inflation model of the black hole 
interior and discuss its limitations. The boundary con- 
ditions will be discussed in section §4. These results 
will be applied to a null crossflow model in section §5. 
In §6, we present an approximate analytic solution to 
the scalar field equations inside a spherical black hole. 

2 Schematic Mass Inflation Mo- 
dels 

The essential physics of mass inflation - a blueshifted 
influx combined with an outflux - can be illustrated 
using simple solvable models. The simplest model is 
that of two concentric spherical null shells, one ingoing, 
the other outgoing, which cross without interaction. 
This situation was examined for shells crossing in the 
Schwarzschild geometry pa and extended to Reissner- 



Nordstrom |16 



Imagine the spacetime split into four regions A, B, 
C and D by the two crossing shells (see Fig. 1 of |1(| ) . 
Each region i, is described by a Reissner- Nordstrom 
solution with metric function fa = 1 — 2?n i /r + e 2 /r 2 . 
The masses of the ingoing and outgoing shells (before 
collision) are mj„ = ra^ — tub and m out — mc ~ 
tub ■ The final result of the analysis is that the metric 
functions are related by the generalised DTR relation 

f A fB=fcfn, (1) 

which holds at ro, the radius of crossing. Mass infla- 
tion occurs when the ingoing shell is very close (within 
distance e) of the Cauchy horizon corresponding to the 
extension of sector B, the space between them before 
the shells cross. Then Jb ~ e, while the product fcJD 
stays finite. This forces the metric function in region 
A, after the shells cross, to diverge as /a ~ 1/e- 

The crucial role of the outflux is illustrated by writ- 
ing the expression for the mass after the crossing, m^ 

H 

mcm D , m m m out . . 

™>a = + — r~r^- <- 2 ) 

m B mB{-jB) 

The first term is finite. When the ingoing shell is close 
to the Cauchy horizon, the second term becomes very 
large. But note that the second term is proportional 



to the product minm ou t. The absence of the outgoing 
shell would render the second term harmless. Thus the 
outflow, however weak, is crucial for the divergence of 
the mass function after the shells cross. 

Ori has introduced a generalization of this simple 
model M. He considers a continuous influx and treats 
the outflux as a thin shell. This allows the match- 
ing of two ingoing Vaidya solutions along the outgoing 
lightlike shell E, a finite Kruskal time after the event 
horizon (see Fig. 1 of p7|): 



ds 2 = dv±(f±dv± - 2dr) + r 2 dtf, 
f± = l-2m ± (v±)/r + e 2 /r 2 , 



(3) 



where the subscript + (-) refers to the region after (be- 
fore) the shell. The Einstein equations link the mass 
with the influx 



L ± = dm±/dv ± , T^ 



± = L ± (v±) 
b 4vrr 2 



d a v±d b v±. (4) 



Continuity of the line element and the radial coordi- 
nate, r, yields the equations 



f+dv + = f-dv- = 2dr 



(5) 



along the shell. Continuity of the influx across the shell 
gives the equation 



1 dm + 1 dm. 



fi dv+ ft dv- 
These two equations yield the simple equation (Ti 
dm + dm^ 

17' 



/- 



(6) 



(7) 



in which mass inflation is evident. The metric function 
/_ goes to zero as the Cauchy horizon is approached, 
causing the right hand side of the equation to diverge. 
The presence of the outgoing shell displaces the ap- 
parent horizon to smaller radii so that /+ ^ at the 
Cauchy horizon. This equation implies that beyond 
the shell, the mass will diverge at CH. 

This model can be solved asymptotically close to 
the Cauchy horizon. The mass function prior to 
the shell must reproduce the Price power law tail: 



m - 



Jl 



3jt(ko^-) ^ p 1%> - Here u_ is the usual 



KO(P-I) 

Eddington-Finkelstein advanced time coordinate which 
is infinite on the Cauchy horizon, f3 is a dimensionless 
constant and kq is the surface gravity of the inner hori- 
zon. Equation (|5|) is then integrated for r along E as 

v_ — ► oo 



rs(u-) 



'•() 







roi4(p- I 
( K oV-)- {p - 1) I I 



(8) 



(P-1) 



Equation (R) can now be integrated using (ph to show 
that the mass function diverges exponentially 



m + {v-)~e KoV -(n V-)- p , 



(9) 



This phenomenon has been dubbed mass inflation J7J. 
Indeed this is a scalar curvature singularity since the 
Weyl curvature invariant diverges exponentially, ^2 ~ 
e K ° v - /tq as the Cauchy horizon is approached. 

It is worth noting that there are coordinates in which 
the metric is finite. Near the CH the line element is 
given to a very good approximation by 



ds 1 



2^±(rdr + m+(v + )dv + ) + r 2 dfl 2 (10) 

r 



where u+ is the standard advanced time coordinate be- 
hind E. It is easily checked, that the coordinate u, 
defined by 

du = rdr + m+(v+)dv+ 

is regular at the CH. (||) now becomes 

ds2 = 2 dv ± du +r2dQ2 

r 

The mass inflation singularity, though much stronger 
than a whimper singularity is still very weak in this 
sense. It is exactly this property of the mass inflation 
singularity, that will allow us to construct an approx- 
imate solution to the full cross flow and scalar field 
equations. 

3 The Continuous Crossflow 
Model 

The previous idealized models show that when the out- 
flux is concentrated as a shell, the mass function (and 
the Weyl curvature invariant) inflates exponentially. 
Originally, this was shown for a continuous, arbitrary 
outflux which starts a finite time after the event hori- 
zon. In this section we shall introduce the notation and 
Einstein equations for spherical symmetry, review the 
standard mass inflation solution and discuss its limita- 
tions. 

In a general spherical spacetimc it is convenient to 
introduce a coordinate system x a — (x a ,9,(/)), (a = 
1,2), where x a are the coordinates of the radial two- 
spaces (8, (f>) = constant. Introducing a function r(x a ) 
which measures the area of the two-spheres, the metric 
can be written: 



ds 2 =g ab dx a dx b + r 2 dtt 2 1 



(11) 



where dft 2 is the line element of the unit two-sphere 
and g ab is the metric of the radial two-spaces. 



by 



Scalar fields m{x a ),f{x a ) and n(x a ) can be defined 



/ = g ab d a rd b r = l- — + -, 
1 



1 e 2 

r,drf = 2"( TO )■ 

2 r z r 



The Einstein field equations are then written as 

G a /3 — 87r(_E Q/3 + T a p), (13) 

where E a p is the Maxwellian contribution to the stress 
tensor. For a point charge of strength e located at the 
origin, this is 



Ei 



87rr' 



-diag(-l, -1,1,1) 



(14) 



(We shall assume that e is fixed.) The non-Maxwellian 
contribution will be decomposed as 

% a = T b a , T*=T+ = P (15) 

where P is the tangential pressure. The energy- 
momentum conservation equation is 

{r 2 T ab ), b = {r 2 y a P. (16) 

The field equations can then be written as 

r-ab + ugab = -~Airr(T a b - gabT) , T — T° 
R-2d r n = 8tt(T-2P). (17) 

Equations (|lj) and ( |17|) yield the equations 

m, a = 47rr 2 {T b a ~8 b a T)r, b (18) 

which can be used to derive a wave equation for the 
mass function 

□ m = 4:n{2rf(P-T) + nr 2 T) - (4tt) 2 r 3 T a6 T ab . (19) 

In order to be more specific, we will write the two- 
metric using null coordinates U, V 



g ab dx a dx b = -2e 2a dUdV. 



(20) 



We will choose V, U to be Kruskal coordinates in which 
the asymptotic Reissner-Nordstrom metric (mass m 
and charge e) is regular on the inner horizon. They are 
related to the usual Eddington-Finkelstein advanced 
and retarded times v, u by 



U 



\ V 



(21) 



where kq = k(tq) is the surface gravity of the Cauchy 
horizon and r$ is its radius. The wave operator acting 
on a scalar ip takes the form 



mp = -2e- 2(y ij tU v. 



(22) 



In this coordinate system, the Ricci scalar is R = 
-2Dcr. 

Null radial vectors pointing inwards and outwards 
can be defined 

l a = -d a V, n a = -d a U. (23) 

In the original mass inflation analysis |7]], a null cross- 
flow stress tensor was used to model the gravitational 
radiation. This effective stress-energy model |Lq] is jus- 
tified by the high blueshift near the Cauchy horizon. 
The stress tensor for null crossflowing radiation can 
then be written as 



T„ 



L in {V) 

47rr 2 



lah 



47rr 2 



-n a i%b 



(24) 



which satisfies the conservation equations (|16[ ) and has 
P = T = 0. The conservation equations force Li n 
(L out ) to be a function only of V (U). 

In the Kruskal coordinate, V, the Price power-law 
tail has the form 



L in (V) 



dm 



P 



dv K dV' ' {-kqV) 2 



(-H-v)Y 



(25) 

As the Cauchy horizon is approached, in the limit, 
V — > 0_, Li n diverges and the source term in the wave 
equation for m diverges as well. The integral solution 
for the mass function is |7J] 

771(17, V) = 

rU r V 

r'- l e- 2 *' L in {V')L ov , t {U')dU'dV' 

'Ui JVi 

+m in (V) + m out (U) - m\. (26) 

The gravitational wave tail influx is turned on at ad- 
vanced time V\ and the outflux is assumed to be 
switched on at the advanced time U\ which is behind 
the event horizon. The divergence of L in (V')dV leads 
to mass inflation with the mass function behaving as 
m ~ 1/V. Of course, this is only true if the combina- 
tion r _1 e~ 2cr does not go to zero, but this was proved 
by Poisson and Israel JTJ]. 

The previous mass inflation analyses suffer some lim- 
itations. In the picture presented Ja p| it is always as- 
sumed that the outflux is turned on abruptly after some 
finite time behind the event horizon. Essentially, this 
amounts to the assumption that a null portion of the 
Cauchy horizon exists because the solution before the 
outflux begins is the Vaidya solution. This segment's 
existence depends on the form of the outflux crossing 
it. If the outflux at early retarded times (U — ► — oo) 
is too strong, a spacelike singularity will form. The 
effect of an outflux crossing a null ray is described by 
Raychaudhuri's equation 



a r 



-47rrT; 



A A 



(27) 



where A is an affine parameter on the null hypersurface 
and T\\ = T ab ^jj-^j- is the transverse flux and dx b /d\ 
is tangent to the generators of the null hypersurface. In 
the case of interest, the null hypersurface is the Cauchy 
horizon and A — > — oo corresponds to its "meeting" 
with the event horizon at H in Fig. 1. 

Examination of ( |27J ) shows that in order for r to be 
finite as A — > — oo , T\\ must satisfy 



X 2 T 



xx 



as A 



(28) 



To test this condition we need to relate the affine pa- 
rameter A to the null coordinate U by 



dU 
~d\ 



JJV 



-2a 



(29) 



If, 



diverges as V — » and U — ► — oo, then depend- 



ing on Tjju, the Cauchy horizon may not survive. 

The earlier models of mass inflation do not address 
this issue since the outflux is turned on after the event 
horizon. The corner region, V — ► 0, — oo < U < U\ 
is described by a Vaidya solution. (The metric func- 
tion e~ 2a is finite in this region for Vaidya.) For more 
general models which include the corner region, the 
behaviour of a must be found. To do this, we need 
to specify the appropriate initial conditions for Tjju 
which arc physically reasonable. 

4 The Outflux 

A star collapsing through its event horizon provides 
two sources of outflux. First, the star shines as it col- 
lapses and will irradiate the Cauchy horizon after the 
event horizon is passed. While we will not attempt to 
describe the actual form of the the star's radiation, we 
do know that in a freely falling frame at the event hori- 
zon, the radiation must be measured to be bounded. 
Kruskal coordinates for the event horizon, [/+ = e K+u 
are appropriate for freely falling observers. These ob- 
servers measure Tu + u + < constant. Transforming this 
to the Kruskal coordinate U appropriate to the inner 
horizon the outflux across CH is 



(Tuu)star = Tu+Uj, 



fdU^ 



J+u+ V dU 



(30) 



As we shall see, the outflux due to the star has a neg- 
ligible effect compared to the backscattering of the in- 
coming radiation. 

Consider the evolution of a massless spherically sym- 
metric scalar field in the black hole interior. The char- 
acteristic initial value problem is completely specified 



by data given on the the event horizon. The physi- 
cal initial data are determined by the Price power law 
wave tail v~ p . 
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Fig 2. View of the spherical black hole interior 
(one angular variable suppressed) , showing fu- 
ture light cones and a stream of infalling radi- 
ation, partially scattered off the potential bar- 
rier, with the remainder accumulating along 
the Cauchy horizon. The point H of Fig. 1. 
now splits into two points Hi and H2. 

For reference purposes, consider the evolution of a 
scalar field in a fixed Reissner-Nordstrom background 
with mass mo- The far right hand side of Fig. 2 (with 
all fluxes turned off), describes the static Reissner- 
Nordstrom solution, ft is distinguished by an outer 
layer where the gravitational potential barrier is weak 
and perturbations can propagate without impediment. 
The potential is peaked around the radius r = e 2 /mo. 
This is where most of the perturbation will be scat- 
tered. Much further in, near r ~ ro, the Cauchy hori- 
zon is approached and infalling radiation is strongly 
blucshiftcd. ft is important to note that the radiation 
is scattered long before it reaches the large blueshift 
zone. This motivates our treatment of the evolution 
of fields as a scattering problem on a static Reissner- 
Nordstrom background. 

Mathematically, scattering of a massless field is given 
by □ "J = 0, which using the usual advanced/retarded 
coordinates is 



V(r) 



V{r)ct> , 
Mr) d 

Ar dr 



vf 



fs(r) 



(31) 



where / has been defined in ( |l2| ) and the subscript s 
denotes the static Reissner-Nordstrom case with mass 



m . 



This is a one-dimensional scattering problem, ft is 
greatly simplified by the fact that the potential V(r) is 
highly localized near e 2 /mo. It falls off exponentially 
H in the tortoise coordinate defined by dx = dr/f s (r) 



near the event horizon x = —00 and the Cauchy hori- 
zon x — 00. A scalar field will propagate freely near the 
event and Cauchy horizons and will only strongly inter- 
act with the curvature in the thin belt around e 2 /mo- 
At the horizons the scalar field solutions will be of the 
form of ingoing and outgoing waves e~ %UJV and e~ tu}u . 
The effect of the potential will be to alter the ampli- 
tudes by the reflection and transmission coefficients, 
r(u) and t(oj). 

If the initial value on the event horizon is 4>o( v ) then 
its Fourier transform Iq] 



4>o(uj) 



I 



>2ir J- 



(j>o{v)e %ulv dv 



(32) 



allows us to write the form of the scattered waves as 
X(v) + Y(u) where 



X(v) 
Y(u) 



<p (uj)t(uj)e- lujv duj 
Mu)r(uj)e luju duj. (33) 



The initial conditions are <j>o{v) = (kov)~ p ^ 2 Q(v — vi) 
where the influx is assumed to start after v\. The 
Fourier transform behaves as |o 



4>o(uj) 



,p/2-i 



(34) 



This can be used to calculate the transmitted and re- 
flected flux. The stationary phase approximation can 
be used to evaluate the integrals (p3). For large v the 
transmitted flux has the form 

X(v) ~ t(LUo)(k v)-P/ 2 , 
u>o = —ip/2v (35) 

t(uj ) ~ l/v 

and the reflected outflux is for large negative u 

Y(u) - rM(-«o«r p/2 , 
wo = —ip/2u (36) 

r(u>o) ~ constant. 

We calculated the reflection and transmission coeffi- 
cients shown in ([35]) and (j^) using a simple model for 
the scattering potential: a rectangular well adjacent 
to a rectangular barrier. Since the actual scattering 
potential looks approximately like a well adjacent to 
a sharp-edged barrier, this model captures the essen- 
tial features of the potential. For low energy scattering 
it is expected that the perturbations will be strongly 
influenced by the potential so that there will be an al- 
most total reflection. Actually, it is more appropriate 
to use the term refraction here since the reflected beam 
continues on to smaller radii. 



The general effect of the Reissner-Nordstrom curva- 
ture is to scatter the influx T vv ~ (kqv)~ p into a re- 
flected flux T uu ~ a (—kqu)~ p and a transmitted flux 
T vv ~ /3(ko^) _p_2 near the Cauchy horizon, where a 
and /3 are the reflection and transmission coeffients and 
are O(l). Is this linearized scattering theory useful for 
our problem? 

Consider the initial layers just beneath the event 
horizon r = r + — e. This region is far above the 
strongly blueshifted region, so the flux of energy is only 
that of the Price gravitational wave tail v~ p . For late 
times this is very weak and will only be a small pertur- 
bation from the vacuum Reissner-Nordstrom solution. 
The effect of the Reissner-Nordstrom geometry on the 
wave tail influx will be negligible and the influx from 
the event horizon will freely propagate to the Reissner- 
Nordstrom potential barrier. 

The black hole interior can be approximated from 
the event horizon to the scattering potential as Reis- 
sner-Nordstrom . Our method will be to solve the Ein- 
stein equations in the interior region after the scatter- 
ing potential. Initial conditions can be set just after 
the potential barrier, given by the Reissner-Nordstrom 
scattering problem. In fact, just after the potential 
barrier the fluxes are not particularly large since the 
blueshift region has not been approached yet. Until the 
Cauchy horizon is reached the perturbations to Reis- 
sner-Nordstrom are small. Our approach is to find an 
approximate solution which is good close to the Cauchy 
horizon (where perturbations are large) and which sat- 
isfies the initial conditions given by scattering from the 
potential barrier. 

As a first step toward the analytic scalar field ap- 
proximation, we shall first introduce a null cross flow 
solution which incorporates the boundary conditions 
discussed here. 



Analytic Approximation 
light like Crossflow 



for 



Our aim is to construct an analytic model for the black 
hole interior after the potential barrier. We shall first 
start with a null crossflow stress tensor and introduce 
some approximations. This can then be used as a 
model for what we expect to happen in the scalar field 
evolution. 

The null crossflow stress tensor is of the form of equa- 
tion (E4f) with the luminosity functions given by scat- 
tering 



L in (V) 

L out {U) 



f3(- Ko V)- 2 (-\n(-V))- q 

a(- KQ uy 2 (H-u))- p 



where a and f3 are dimensionless positive numbers, cor- 
responding to the reflection and transmission coeffi- 
cients respectively and q = p + 2. 

Introduce the functions A(U),B(V) defined by 



L out {U) = A"{U) 
L in {V) = B"{V) 



(38) 



where ' denotes ordinary differentiation. For V — > 



B(V) = 

B'(V) = 

and for U — ► — oo 
A(U) = 

A\U) = 







«§(«-!) 



(-H-v)) 



-9+1 



1 







«§(-V) 



(-ln(-V))- 



-H-v) 



(39) 



«8(P-1) 



(ln(-L0)- p+1 



1 



H-u) 



K*(-U) 



1 - 



(H-u))- 



H-u) 



(40) 



In the corner region V — ► 0, U — > — oo, the functions A 
and B are small, but derivatives of B with respect to 
V diverge. 

We wish to concentrate on the region after the poten- 
tial barrier at early times. Before the potential barrier 
we expect Reissner-Nordstrom to be a good model. In 
the innermost regions we must model the effect of the 
infinite blueshift of the inflowing radiation. 



Using the metric (20 ) we note that the Einstein equa- 
tions allow us to write wave equations for two combi- 
nations of the metric functions, which do not include 
the mass function as a source term: 



(\n(r- l e- 2 °)),uv 
{r 2 ),uv 



e 2a 
2^2 

„2cr 



1 



1- 



3V 

r ;2 



(41) 



(37) 



In order to solve the evolutionary problem, we need 
to find the solution in the intermediate region after the 
potential barrier and before the region where r goes 
to zero. This region will be defined by r ^ 0. When 
this stipulation is made it is impossible for r~ 1 e -2cr to 
go to zero |0. This means that both wave equations 
mWj do not have any potentially diverging source terms 



and both will have finite solutions. For conciseness, 
introduce the bounded and non-zero variables 



the Einstein equations to find the first order correction 
terms. Equations (E3J) can be integrated to solve for p: 



X = r~ l e- 2 ° , p 



1 



(42) 



P = Ps 



(50) 



The Einstein equations can then be written as equa- 
tions (Eh) and the null hypersurface constraint equa- 
tions: 



dV" 



v" 



X (U",V") 

U du „ r U" 



X {U",V") 



dV'x(U',V')B"(V') 
dU' X (U',V')A"(U'). 



du(xp,u) 

dv(xP.y) 



= -XA" 
= -XB". 



The mass function obeys the wave equation 
m,vv = xA"B". 



(43) 



(44) 



As long as r ^ 0, we can write a solution with x 
and p being close to their Reissner-Nordstrom values 
plus perturbations which are small in this region. The 
metric functions for static Reissner-Nordstrom with a 
mass mo will be denoted with a subscript "s" , so that 
f s (r s ) and n s (r s ) are defined by equation (|12|). The 
functions p s and Xs and their derivatives are 

1 2 

Ps = ~ 2 r s 



Ps,U 



Ps,V 



lr s f s 



2k U 


lr s f s 


2k V 


-2UV kl 



(45) 



It is clear in our approximation scheme that ( |50| ) is 
the leading order contribution to the solution of the 
Einstein equations. The contribution from ( |4l| ) will be 
of lower order. 

Integration of (|50|) by parts gives the solution 



p = Ps - {A + B) + e, 



(51) 



where e is 



dV" 



Xs{U",V") 



v" 



dV' Xs y'(U',V')B'(V') 



dU" 



u" 



dV'xs,u'(U'X)A>(U>) 



Xs{U",V") 

V I BdV + V f AdU (52) 



which is much smaller than A + B in the remote past of 
CH. Using this approximation in the second equation 
of <MW) and expanding to first order in A and B, allows 
the estimation 



X = Xs + 0{^{A + B)). 



(53) 



Xs,U 

Xs,v -- 



k V ( 2n r s 

I -I : (k - K a ) 



' s 
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Is 

2n Q r s 
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(«0 - «a) 



/. 



1 



-Jo, Pa,U 



-2UV 

r V 



«•■() 



", Ps,V 



X* 



I'Q 



Xs,U 



K V 

' — — ' ^y 



r U 

K U 



(46) 



In the limit of the Cauchy horizon, (UV — > 0, r s — > r ) 
these functions take on the limiting value 



(47) 
(48) 

(49) 



We can now construct a solution to the Einstein 
equations using an iterative approach, taking the static 
Reissner-Nordstrom solution as the zeroth order solu- 
tion (x' ' = Xs,P = Ps) and substituting back into 



Substitution of this order of correction back into ( p0[ ) 
yields a second order approximation 

P = p s - A(l + 0(ln(-U))- p ) 

- B(l + 0{-\n{-V)y q ). (54) 

Clearly, in the corner region where (hi(— U))~ l ~ 
(— ln(— V))~ ~ 0, p is well approximated by the lead- 
ing order solution (pl[). 

To linear order in A and B, the mass function can 
be integrated from ( J44| ) 

m(U, V) = XsA'B'{\ + 0(A + B)) + m' m + m' out - m 

(55) 
which in the limit V —> is 

m ~ ^2^rAH-u))- p {-H-v))- (l (56) 

r /vg UV 
showing the usual 1/V inflation found in earlier work 



The solutions for the original metric functions r and 
a are 

(A + B) 2AB „, i9 9N 

r = r, -^ '- + — s- + 0{A 2 + B 2 ) 



a = a. 



r s r% 

(A + B) AB 
2r? + 7f 



0(yl 2 + S 2 ). (57) 



This approximation is not so good as the scatter- 
ing surface is approached (UV — > 1, so that correction 
terms (p2) are comparable to the first order terms in 
d5j])). We already know that the solution near the 
scattering surface should be approximately described 
by the Reissner-Nordstrom solution. It is only after 
this region, deep into the blueshift region that an ap- 
proximate solution is needed and this is where it is 
important that the solution be accurate. The solution 
that we have found is accurate where it matters, close 
to the Cauchy horizon. 

6 The Scalar Field Solution 

Using approximations similar to those just discussed 
for lightlike crossflow, we can develop an approximate 
analytic solution for the scalar field equations. As be- 
fore the physics tells us that the interior solution can 
be approximated well by the static Reissner-Nordstrom 
solution from the event horizon, down until the scat- 
tering surface. 

The Einstein equations for coupling to a massless 
scalar field are 



r<p,uv 
P,uv 



(Mx)). 



UV 



-r,u<P,v - t,v<P,u 

r\ V r 
in(j)u4>,v 



(58) 



1 / 3£ 

r 3 x \ r 2 



(xp,u),u = -^px4> 2 u 

{XP,v),V = Snpxpy 

m : uv = ^n 2 xr i <f?u<t?v-^ r f ( t>,u(i>,v- 

Define functions a(U) and b(V) by setting their 
derivatives equal to 



a'(U) = <f>,u\ h , V(V) = <t>,v\ b 



(59) 



where the subscript b refers to the value of the scalar 
field given by scattering at the underside of the poten- 
tial barrier. 

As before, the wave equations for p and \ have solu- 
tions which are finite and non-zero in the corner region 
as long as <j) does not diverge. The initial conditions 
given by scattering (Ba,Bq) are that the scalar field is 



initally regular. Near the scattering surface the radius 
will be close to its Reissner-Nordstrom value, so using 
the scalar wave equation and the Reissner-Nordstrom 
radius (|4q), it can be seen that the UV mixed deriva- 
tive of the scalar field, near the initial surface is 



;uv\b 



1 /. 



?,v 
U 



<t>,i 



01* 



2r s K U V 

- (- ln(-V)y q/2 + (ln(-[/))- p / 2 . (60) 

This derivative is small in the corner (V — > 0, U — > 
— oo), so in the earliest regions the scalar field will not 
be changing rapidly from its initial value. This mo- 
tivates us to make the ansatz that the leading order 
behaviour of the scalar field, near the Cauchy horizon 
should be 



$ 



^ = b' 



I A" 
Airr 2 



-TT\\-p/ 2 



I B" 



^rlnli-V) 



(H-u)) 



(61) 



(-H-v))-^ 2 



With this ansatz, we can see that the scalar field will be 
small everywhere in the corner region, but that deriva- 
tives with respect to V will diverge near the Cauchy 
horizon. 

As before we can calculate the first order correction 
terms by iterating the Einstein equations, again taking 
the zeroth order solutions for p and x to be the same as 
the Reissner-Nordstrom solution. The solution for p is 
the same as the lightlike cross flow solution (51). Sub- 
stitution of ( |5l| ) and (|6l]) into the scalar wave equation 
yields the first order equation 

<f>,uv - -^-((Ps.u - A')b' + {p sy - B')a') (62) 

2/3 s 

which can be integrated asymptotically in the corner 
region, making use of rt45j) 



-*(Ab + Ba) + OCUVia 



Using the first order solutions for 
equation for x can be integrated 



6)). (63) 
and p the wave 



In x = In Xs + 8irab+0(A + B). (64) 

To leading order the mass function is integrated to be 

m(U, V) - kIjtqA'B' (65) 



and the metric functions r and a are 
(A + B) AB 



0{A 2 + B A 



cr. — Anab - 



0(A 2 +B 2 ) 



{A + B) AB 

2r 2 + ~rj 



(66) 



The existence of the Cauchy horizon in this solution 
can now be tested. Substitution of the solution for a 
given by ( |66| ) into (p9|) gives the following asymptotic 
relation between the affine parameter, A, and the co- 
ordinate, U, 



A = -^(1 + 0(a)) ,U 



(67) 



Condition (g8|) for the scalar field solution now reads 



lim A 2 T A a = 



(68) 



lim 



M-u)) 



-p+2 



Anr^ 



(1 + 0{ln(-U))- p+1 ) =0. 

Since condition (|2g| ) is satisfied, the Cauchy horizon ex- 
ists in our approximate solution to the Einstein-scalar 
field equations. This is of course evident directly from 
the asymptotic form of the metric given by (|2C| ) and 

7 Conclusions 

We have found an approximate solution to the scalar 
field equations in the black hole interior which includes 
a complete past segment of the Cauchy horizon and 
contains the requisite number of arbitrary functions to 
be considered "general". In our picture (see Fig. 2) 
the Cauchy horizon will be focussed by the outgoing 
flux, forcing it to eventually contract to r = 0, forming 
a singularity which is possibly spacelike. But we find 
that the r = singularity does not precede the Cauchy 
horizon in the earliest region of the interior. 

This is not in agreement with a previous numerical 
study pi| . That work consisted of two simulations, in 
the first of which |H| initial conditions were placed on 
the event horizon. The result of this simulation was 
that the Cauchy horizon survived the introduction of 
a scalar field. In the second study, initial conditions 
were set outside the event horizon, which in effect al- 
lows tails to form. These results (see their Fig. 4) 
show that any null portion is abbreviated or perhaps 
completely absent. It is difficult to judge from their ev- 
idence, but the authors state that the Cauchy horizon 



is destroyed and replaced by a spacelike r = singu- 
larity. Further numerical work on this problem is in 
progress [|20| , which should help clear up the discrep- 
ancy in results. 
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